Model  Development  and  Model-Based  Control  Design 
for  High  Performance  Nonlinear  Smart  Systems 

AFOSR  FA9550-04-1-0203 


Final  Report 

for  the  period 

April  1,  2004  -  September  30,  2007 


Program  Manager: 

Scott  Wells 

AFOSR/NM 

4015  Wilson  Blvd.,  Room  713 
Arlington,  VA  22203-8409 
scott.wells@afosr.af.mil 
703-696-7796 


Principal  Investigator: 

Ralph  C.  Smith 

Center  for  Research  in  Scientific  Computation 
Department  of  Mathematics 
North  Carolina  State  University 
Raleigh,  NC  27695-8205 
rsmith@eos.ncsu.edu 
919-515-7552 


i 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

20-11-2007  Final  Report  01SEP2004-30SEP2007 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 

FA9550-04-1-0203 

Model  Development  and  Model-based  control  design  for  high  5b.  GRANT  NUMBER 

performance  nonlinear  smart  systems 

5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 


North  Carolina  State 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFOSR 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Distribution  A:  Approved  for  Public  Release 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL-SR-AR-TR-08-0060 


14.  ABSTRACT 

We  have  developed  a  unified  energy-based  framework  for  quantifying  hysteresis  and 
constitutive  nonlinearities  inherent  to  piezoelectric,  magnetic,  and  shape  memory  compounds 
which  is  amenable  to  inversion  and  subsequent  use  as  an  inverse  filter  for  linear  control 
designs.  We  have  also  developed  stochastic  modeling  techniques  which  quantify  the  highly 
complex  stiffness  properties  of  ionic  polymers  in  a  manner  which  facilitates  design  in 
applications  ranging  from  biological/chemical  detection  to  robotic  designfor  aerospace 
structures.  The  control  component  has  focused  on  the  development  of  robust  linear  designs 
exploiting  nonlinear  filters  and  fully  nonlinear  algorithms  which  incorporate  modeled  physics 
directly  into  the  control  design. 

Open  loop  control  experiments  have  been  performed  and  present  investigations  are  focused  on 
closed  loop  experimental  validation  of  the  control  theories. 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


Objectives 

This  research  program  focused  on  the  development  of  energy-based  and  stochastic  models,  numerical 
approximation  techniques,  and  model-based  control  designs  for  piezoelectric,  magnetic,  ionic  poly¬ 
mer,  and  shape  memory  alloy  actuators  and  sensors  operating  in  highly  nonlinear  and  hysteretic 
regimes.  The  first  component  of  the  program  focused  on  the  characterization  of  hysteresis,  consti¬ 
tutive  nonlinearities,  and  thermal,  stress  and  frequency-dependencies  in  a  manner  that  promotes 
transducer  and  control  design.  Secondly,  we  developed  linear  and  nonlinear  control  techniques  that 
exploit  the  physics  encapsulated  in  models  to  address  stringent  control  objectives  (e.g.,  micron-level 
tracking,  kilohertz  to  megahertz  operating  regimes)  while  operating  in  highly  nonlinear  regimes. 

Accomplishments 

Actuators  and  sensors  comprised  of  piezoelectric  (PZT),  magnetic,  shape  memory  alloy  (SMA),  and 
ionic  polymer  components  provide  unique  design  and  control  capabilities  for  a  range  of  aeronautic 
and  aerospace  applications.  Piezoelectric  compounds  are  lightweight  and  provide  broadband  sen¬ 
sor  and  actuator  capabilities.  Furthermore,  geometry  enhancement  of  outputs  in  architectures  such 
as  THUNDER  (THin  layered  UNimorph  Driver  and  sEnsoR)  transducers  provide  large  strain  or 
displacement  capabilities.  Due  to  these  attributes,  PZT-based  devices  are  presently  considered  for 
synthetic  jets,  shape  modification  of  airfoils,  high-speed  valve  and  fuel  injector  design,  and  energy 
harvesting.  SMA  devices  provide  the  highest  work  densities  of  the  previously-mentioned  materials 
but  are  presently  limited  to  frequencies  below  100-200  Hz.  Within  the  context  of  Air  Force  appli¬ 
cations,  SMA  are  being  considered  for  chevron  design  to  reduce  jet  noise  and  increase  performance, 
deployment  and  vibration  attenuation  in  modular  antennas,  and  shape  modification  of  airfoils  [22]. 
Magnetic  actuators  provide  large  force  capabilities  for  applications  including  remote  structural  de¬ 
ployment  as  well  as  high  fidelity  sensing.  Finally,  ionic  polymers  are  under  development  for  use  in 
applications  ranging  from  biological/chemical  sensing  to  remote  robotics.  However,  all  of  these  com¬ 
pounds  exhibit  highly  nonlinear  dynamics  which  must  be  incorporated  in  models  and  model-based 
control  designs  before  their  potential  can  be  fully  realized  in  Air  Force  applications. 

Development  of  a  Unified  Modeling  Framework  for  Ferroic  Compounds 

A  hallmark  of  this  program  was  the  development  of  a  unified  characterization  framework  for  fer¬ 
roelectric  (e.g.,  PZT),  ferromagnetic,  and  ferroelastic  (e.g.,  SMA)  materials  —  which  are  collectively 
designated  as  ferroic  compounds.  As  detailed  in  [22,29],  the  framework  is  constructed  in  two  steps. 
In  the  first,  Helmholtz  and  Gibbs  energy  relations  are  constructed  at  the  mesoscopic,  or  lattice,  scale 
to  quantify  the  local  average  polarization,  magnetization,  or  strain.  In  the  second  step,  stochastic 
homogenization  techniques  are  invoked  to  construct  macroscopic  models  that  incorporate  the  ef¬ 
fects  of  material  nonhomogeneities,  polycrystallinity,  and  nonuniform  effective  fields.  The  resulting 
models  are  sufficiently  accurate  for  material  characterization  and  sufficiently  efficient  to  facilitate 
model-based  control  design. 

We  summarize  here  primary  steps  in  the  framework  and  note  that  details  regarding  the  frame¬ 
work  can  be  found  in  [22-24,27,29].  Extensions  the  theory  to  incorporate  stress  and  temperature- 
dependence  are  provided  in  [1,2,4,6,10,16,21,26,28]. 

Let  P,  M  and  e  respectively  denote  polarization,  magnetization  and  strains  and  let  E,  H  and  o 
denote  an  applied  electric  field,  magnetic  field,  and  stress.  One  can  consolidate  notation  by  letting 
e  =  P,  M,  e  denote  an  order  parameter  and  let  <p  =  E,H,o  denote  conjugate  fields.  Additionally,  we 
let  q  designate  electromechanical  or  magnetomechanical  coupling  coefficients. 
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In  the  first  step  of  the  development,  we  consider  the  Gibbs  energy  relations 


G(E ,  P,  T)  =  ip(P,  T)  —  EP  —  as  Ferroelectric 

G(H ,  M,  T)  =  ip(M,  T )  —  HM  —  as  Ferromagnetic  (1) 

G(ct,  e,  T)  =  ip£(s,  T )  —  ere  Ferroelastic 


where  T  is  temperature  in  degrees  Kelvin.  For  ferroelectric  (e  =  P )  and  ferromagnetic  (e  =  M ) 
materials,  the  the  Helmholtz  energy  relations  have  the  form 


V<e,  T)  =  ipe(e,  T)  +  ^Yes2  -  qs(e  -  eR )2 


(2) 


where  the  piecewise  quadratic  component  ipe  incorporates  the  bistable  behavior  associated  with 
equilibrium  dipole  or  moment  configurations  as  depicted  in  Figure  1(a).  The  second  and  third  terms 
on  the  right  side  of  (2)  incorporate  the  elastic  and  electromechanical  or  magnetomechanical  coupling 
energy.  For  ferroelastic  materials,  V’e  characterizes  the  tristable  martensite-austenite  phase  behavior 
shown  in  Figure  1(b). 

For  regimes  in  which  thermal  activation  is  negligible,  piecewise  linear  kernel  relations  of  the  form 


p=  £  +2<;£-A,'pa. 

—  2  qs  r]  —  2  qs 


(3) 


where  A  =  ±1,  are  obtained  through  enforcement  of  the  equilibrium  conditions  =  0.  The  relations 
for  M  and  s  are  similar. 

To  incorporate  thermal  activation,  which  is  manifested  at  the  macroscopic  level  as  creep,  relax¬ 
ation,  reptation,  or  accommodation,  it  is  necessary  to  balance  the  Gibbs  energy  with  the  relative 
thermal  energy  kT /V  via  the  Boltzmann  relation 

H(G)  =  Ce~GV/kT.  (4) 


Here  k,  V  and  C  respectively  denote  Boltzmann’s  constant,  a  reference  volume,  and  an  integration 
constant.  This  yields  kernel  or  hysteron  relations  of  the  form 

e  =  x+  (e+)  +  x -  (e_) 

(5) 

s  =  x+  (e+)  +  x -  (e_)  +  xA  (eA) 

where  x+,X-,xA  are  moment,  dipole  or  phase  fractions  and  the  bracketed  quantities  denote  av¬ 
erage  polarizations,  magnetizations,  or  strains  associated  with  specified  dipole,  moment  or  phase 
configurations. 

In  the  second  step  of  the  development,  the  effects  of  polycrystallinity,  material  nonhomogeneities, 
and  variable  conjugate  fields  are  incorporated  by  assuming  that  certain  material  properties  (e.g., 


Figure  1:  Helmholtz  energy  for  (a)  ferroelectric  or  ferromagnetic  materials,  and  (b)  ferroelastic 
materials. 
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coercive  fields  EC,HC,  interaction  fields  Ej,P[eoi  and  relative  stresses  aR)  are  manifestations  of 
underlying  distributions  rather  than  constants.  We  denote  corresponding  densities  by  v\  and  V2- 
Enforcement  of  the  equilibrium  conditions  ^  =  0  and  incorporation  of  Kelvin-Voigt  damping 
subsequently  yields  the  constitutive  relations 

!roo  ro o  _ 

[P(E,T)\(t)=  /  ^{E  +  E^T-E^E^E^dE^Ec 

Jo  J-0 o  (6) 

a  =  Ye  Ed-  q(P  -  PR)2 , 

!roo  roo  _ 

[M(H,T)\(t)=  /  [M(H  +  Hj,T-,  Hc)i/i(Hc)i/2(HI)dHIdHc 

Jo  J-oo  (7) 

a  =  Ye  +  d  —  q(M  —  MR)2, 

{roo  roo 

[e(a,T)](t)=  /  [e((7 +  ai,T-,crR)i/i(crR)v2(cri)dcrRdaR  (8) 

Jo  J-oo 

for  ferroelectric,  ferromagnetic  and  ferroelastic  materials. 

Temperature  changes  due  to  convection,  conduction,  dipole/moment  switching  or  phase  changes, 
and  Joule  heating  can  be  quantified  by  the  general  relation 

cMT(t)  =  -n[hc  +  \/£}[T{t)  -  TE(t)\  +  J(t )  -  ^  haxa.  (9) 

a 

Here  c,  J\4,  hc,fl,  X,£  and  TE(t)  respectively  denote  the  specific  heat  for  the  actuator  material,  the 
mass  of  the  actuator,  a  convection  coefficient,  the  surface  area  of  the  material,  the  thermal  conductiv¬ 
ity  of  the  surrounding  medium,  the  interval  over  which  conduction  occurs,  and  the  time  varying  tem¬ 
perature  of  the  adjacent  environment.  The  components  —hcQ[T(t)  —  TE(t)\  and  —A f2[T(t)  —  TE(t)}/£ 
respectively  quantify  the  effects  of  convection  and  conduction,  and  J(t)  characterizes  Joule  heating 
mechanisms.  The  final  component  quantifies  heat  transduction  due  to  dipole/moment  switching  or 
phase  transitions  so  that  ha  represents  specific  enthalpies. 

In  combination,  relations  (6)  -  (9)  characterize  the  1-D  constitutive  behavior  of  PZT,  magnetic 
materials,  and  SMA,  and  provide  a  framework  for  constructing  distributed  actuator  and  sensor 
models.  The  performance  of  this  modeling  framework  is  illustrated  for  nickel,  PZT,  and  thin  film 
SMA  in  Figure  2  through  a  comparison  and  prediction  of  experimental  data.  This  illustrates  both 
the  accuracy  of  the  methodology  for  a  variety  of  regimes  including  biased  minor  loops  which  ex¬ 
hibit  nonclosure  due  to  thermal  relaxation  or  creep  mechanisms  —  e.g.,  the  nickel  data  and  fit  in 
Figure  2(a)  and  (b). 


Structural  Models 


Electromechanical  and  magnetomechanical  constitutive  relations  of  the  form  (6)  and  (7)  provide 
the  basis  for  constructing  structural  models  for  systems  employing  ferroelectric  or  magnetic  sensors. 
To  illustrate,  consider  a  piezoelectric  (ferroelectric)  rod  having  cross-sectional  area  A,  length  £, 
density  p,  Young’s  modulus  Y,  Kelvin-Voigt  damping  coefficient  c,  a  fixed  end  at  x  =  0  and  a  free 
end  at  x  =  £.  As  detailed  in  Chapter  7  of  [22],  force  balancing  yields  the  relation 


dx 


(10) 


where  u  is  specified  in  (6)  and  u  denotes  the  longitudinal  displacement.  Use  of  the  linear  strain- 
displacement  relation  e  =  and  integration  by  parts  yields  the  weak  model  formulation 


re  .d2u  f e  dud<t> 

pA—=(j)dx+  /  Y  A——dx 
3  OX~  Jo  ox  dx 


+  fCApkfd*=A^p-prff 

J o  oxot dx  J  o 


d<f> 

dx 


dx 


(11) 
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Figure  2:  (a)  Nickel  200  H-M  data,  and  (b)  H-M  model  fit  for  multiple  biased  minor  loops  from  [24]. 
(c)  Model  comparison  to  quasistatic  PZT5H  data  from  [27],  and  (d)  model  fit  to  thin  film  SMA  data 
from  [10]. 


for  all  test  functions  (j)  £  Ffd(0,^)  and  P  given  by  (6).  This  provides  a  framework  appropriate  for 
approximation  and  subsequent  model-based  control  design.  Similar  analysis  is  provided  for  beams, 
plates  and  shells  in  [22]. 

Stochastic  Models  for  Ionic  Polymers 

Ionic  polymers  exhibit  a  highly  complex  structure  comprised  of  a  hydrophobic  backbone  and 
hydrophilic  side  chains  which  presently  precludes  energy  analysis  in  a  manner  that  permits  real-time 
implementation.  During  this  program,  we  have  developed  stochastic  models  at  the  molecular  level 
which  can  be  employed  in  low-order  macroscopic  models  to  predict  effective  material  properties  such 
as  stiffness  [11,30].  By  employing  analytic  density  estimation  techniques,  the  macroscopic  models  are 
sufficiently  efficient  to  permit  both  actuator/sensor  characterization  and  to  consider  the  significantly 
more  difficult  inverse  problem  of  material  design  to  achieve  specified  control  criteria.  This  provides 
an  initial  step  toward  the  development  of  a  paradigm  to  custom  design  materials  to  optimally  achieve 
control  criteria  for  Air  Force  Applications. 

Parameter  Estimation 

A  crucial  step  for  model  construction  or  updating  is  the  development  of  highly  robust  parameter 
estimation  techniques  to  identify  the  densities,  and  other  model  parameters,  based  on  attributes  of 
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the  data  with  little  input  required  from  users.  To  achieve  this,  we  have  developed  theory  and  Graph¬ 
ical  User  Interface  (GUI)-based  packages  that  provide  initial  parameter  estimates  based  solely  on 
measured  properties  of  the  data  [25] .  Representative  fits  to  frequency-dependent  data  are  illustrated 
in  Figure  3.  We  note  that  these  fits  required  approximately  4  minutes  of  CPU  time  on  a  personal 
computer  and  hence  the  algorithms  are  sufficiently  efficient  to  permit  rapid  model  construction. 


-1000  0  1000  2000  3000  4000 

Current  (mA) 


x  10~5 


(a) 


(b) 


Figure  3:  Model  fit  to  magnetic  displacement  data  collected  at  (a)  100  Hz  and  (b)  500  Hz. 


High  Speed  FPGA  Implementation 

High  speed,  high  accuracy  Air  Force  applications  dictate  that  model-based  control  techniques 
be  implemented  at  kilohertz  to  megahertz  rates.  For  the  homogenized  energy  model,  this  means 
that  quadratures  associated  with  the  formulations  (6)  -  (8)  and  construction  of  the  hysterons  (5) 
must  be  performed  in  microseconds.  To  achieve  these  rates,  we  initiated  investigations  regarding 
the  use  of  field  programmable  gate  arrays  (FPGA).  These  user-configurable  hardware  chips  combine 
tens  to  hundreds  of  thousands  of  look-up  tables,  flip-flops,  and  switches  that  can  be  programmed 
by  a  user  to  achieve  a  custom  design.  These  devices  are  naturally  parallelizable  which  proves  highly 
advantageous  when  constructing  hysterons  and  numerically  approximating  the  integrals. 

As  detailed  in  [3,5],  the  use  of  FPGA  yields  speedups  of  approximately  100-fold  for  the  thermal 
relaxation  model  as  compared  with  optimized  C  implementations.  Depending  on  the  number  of 
quadrature  points,  this  yields  implementation  rates  in  the  kilohertz  range  for  inverse  models  which 
is  sufficiently  fast  for  real-time  control  implementation  of  high  speed  Air  Force  applications. 

Model-Based  Control  Design:  Nonlinear  Inverse  Compensators 

The  goal  of  the  control  component  of  the  investigation  was  to  develop  theory  and  algorithms 
that  permit  real-time  tracking  and  vibration  attenuation  using  smart  material  actuators  operating 
in  highly  nonlinear  and  hysteretic  regimes.  For  the  applications  under  consideration,  this  can  result 
in  micron-level  (or  smaller)  tolerances  and  quasistatic  up  to  kilohertz  operating  frequencies. 

As  depicted  in  Figure  4,  there  are  essentially  two  control  strategies  for  hysteretic  systems  using 
nonlinear  models.  The  first  is  to  use  the  characterization  framework  to  construct  an  approximate 
nonlinear  inverse  that  linearizes  the  actuator  response  in  the  manner  depicted  in  Figure  4(a).  Linear 
control  algorithms  are  then  used  to  achieve  control  objectives.  The  second  strategy  is  to  construct 
nonlinear  control  designs  which  yield  input  signals  that  directly  incorporate  actuator  nonlinearities 
as  depicted  in  Figure  4(b).  The  latter  strategy  is  advantageous  when  tracking  trajectories  are  known 
a  priori  whereas  the  inverse  compensators  will  be  required  when  attenuating  unmodeled  or  stochastic 
uncertainties. 
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Nonlinear  Nonlinear 
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(a)  (b) 

Figure  4:  (a)  Linear  control  design  employing  a  nonlinear,  model-based  inverse  compensator  and 
(b)  nonlinear  model-based  control  design. 


During  this  program,  we  developed  highly  efficient  inverse  model  formulations,  based  on  the 
homogenized  energy  framework,  that  provide  a  significant  step  toward  real-time  implementation 
capabilities.  As  detailed  in  [4],  these  inverse  algorithms  employ  either  lookup  tables  or  rational 
function  approximations  in  the  forward  models  in  combination  with  fixed-point  algorithms  to  con¬ 
struct  highly  efficient  inverse  representations.  The  accuracy  of  these  inverse  models  is  illustrated  in 
Figure  5.  Present  C  implementations  of  these  inverse  models  run  at  500  to  1000  Hz  whereas  FPGA 
implementations  are  projected  to  be  approximately  2  orders  of  magnitude  faster. 
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Figure  5:  (a)  Reference  polarization  P,  and  (b)  electric  field  E  given  by  the  inverse  model  with 
parameters  estimated  for  PLZT.  (c)  Comparison  of  the  polarization  P  given  by  (6)  with  input  E 
and  the  reference  polarization  P,  and  (d)  the  absolute  error  \P  —  P\. 
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Figure  6:  Tracking  performance  in  open  loop  experiments  utilizing  a  model-based  inverse  and  in  the 
absence  of  a  filter. 


We  have  also  investigated  numerical  algorithms  and  experimental  implementation  of  open  loop 
designs  employing  inverse  filters  based  on  the  homogenized  energy  framework.  These  results  are 
reported  in  [8,  9,  12]  with  representative  experimental  tracking  results  for  filtered  and  unfiltered 
designs  shown  in  Figure  6. 


Model-Based  Control  Design:  Nonlinear  Optimal  Control 

The  second  main  thrust  of  the  control  component  of  the  program  has  focused  on  the  formulation, 
numerical  implementation,  and  experimental  implementation  of  nonlinear  perturbation  control  de¬ 
signs  for  systems  employing  PZT  and  magnetic  actuators  operating  in  hysteretic  regimes.  As  detailed 
in  [13-15, 17-20],  these  designs  are  developed  in  two  steps:  (i)  construction  of  an  optimal  nonlin¬ 
ear  open  loop  control  signal,  and  (ii)  construction  of  a  perturbation  feedback  component  obtained 
through  linearization  about  the  optimality  system. 

To  illustrate,  let  r  denote  a  reference  signal  to  be  tracked,  y(t)  =  Cx(t )  denote  observations  and 
let  e(t)  =  Cx(t )  —  r(t)  designate  the  error.  The  augmented  penalty  functional  is  taken  to  be 


J  =  \[Cx(tf)  -  r(tf)]TP[Cx{tf)  -  r(tf )]  + 


r*/ 


'to 


7i  —  XT  (t)x(t) 


dt 


(12) 


where  A  denotes  the  adjoint  variable,  the  Hamiltonian  is 


H  =  ^[e{t)T Qe(t)  +  uT(t)Ru(t)]  +  XT  [ Ax(t )  +  [B(u)](t)]  , 


(13) 


and  Q,  R  respectively  penalize  large  errors  and  control  inputs. 

Enforcement  of  necessary  conditions  to  minimize  (12)  yields  the  control  input  relation 


along  with  the  two-point  boundary  value  problem 


(14) 


where  z  =  [x,  A]T  and 


F(t,z) 


z{t)  =  F(t,z ) 


Ax(t)  +  [5(w)](i) 

-ATX{t)  -  CTQCx(t )  +  CTQr(t) 


(15) 


(16) 
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To  approximate  the  solution  to  (15),  we  employ  a  finite  difference  discretization  defined  on  the 
grid  tj  =  jAt,  where  At  =  ^  and  j  =  0,  -  •  • ,  N.  Letting  Zj  ~  z(tj),  this  yields  the  discrete  system 

^  [zj+i  ~  zj\  =  ^  [ F(tj ,  Zj)  +  F(tj+i,zj+1)] 

E0z0  =  [xQ,  0]t  (17) 

EfzN  =  [0,-CTPr(tf)]T. 

The  solution  of  (17)  can  be  expressed  as  the  problem  of  finding  Zh  =  [zo,  •  •  • ,  zn ]  which  solves 

H*h)  =  0.  (18) 


A  quasi-Newton  iteration  of  the  form 


where  solves 


~k+l  _  A  |  tfe 
Zh  —  zh  T  Shi 

E\zkh)eh  = 


is  then  used  to  approximate  the  solution  to  (18).  The  Jacobian  has  the  form 


(19) 

(20) 


So  R0 
Si  Ri 


E\zh)  = 


E0 


Sat-1  Rn—1 


E 


f  J 


(21) 


where 


1 

'  I  o' 

1 

A  ^B[u*}  - 

At 

0  / 

~  2 

-CtQC  -At 

(22) 


The  representation  for  Ri  is  similar. 

It  is  shown  in  [18, 19]  that  an  analytic  LU  decomposition  can  be  determined  for  E(z^).  This 
significantly  reduces  memory  requirements  and  is  fundamental  for  efficient  solution. 

The  relation  (14)  provides  an  open  loop  control  signal  that  is  optimal  for  given  Q,R  in  the 
absence  of  model  or  measurement  error.  To  provide  robustness  with  regard  to  such  uncertainties,  we 
consider  perturbation  feedback.  Linearization  about  the  optimality  system  yields 


6x(t)  =  A6x(t)  +  B5u(t) 
Sy(t )  =  C8x(t ) 


(23) 


where  5u,5x  and  5y  are  first-order  variations  about  u*,x*  and  y*. 

To  facilitate  experimental  implementation,  we  use  classical  PI  control  to  compute  5u;  that  is,  we 
take 

Su(t)  =  —Kpe(t)  —  Kj  [  e(s)ds.  (24) 

Jo 

The  final  control  input  is  then 

u{t)  =  u*{t)  +  5u(t).  (25) 

The  performance  of  this  perturbation  control  in  tracking  experiments,  conducted  at  1  kHz,  is 
illustrated  in  Figure  7.  The  results  in  Figure  7(a)  illustrate  that  an  optimized  PI  design,  in  the  absence 


of  the  nonlinear  open  loop  signal,  exhibits  significant  errors  in  both  the  phase  and  amplitude  due  to 
the  frequency-dependent  hysteresis  shown  in  Figure  7(b).  The  perturbation  control  accommodates 
this  nonlinear  behavior  and  achieves  accurate  tracking.  This  is  an  approximately  30-fold  improvement 
in  operating  frequency  over  previous  control  techniques  based  on  Preisach  models. 

A  crucial  aspect  of  this  approach  is  the  fact  that  the  optimal  open  loop  control  signals  are 
constructed  offline  and  hence  implementation  speeds  are  dictated  by  the  efficiency  of  the  perturbation 
feedback  algorithm.  The  experimental  results  at  1  kHz  demonstrate  the  potential  of  the  technique 
for  applications  that  require  high-speed  tracking  or  broadband  vibration  attenuation  using  actuators 
that  are  operating  in  highly  hysteretic  and  nonlinear  regimes. 


Figure  7:  (a)  Reference  trajectory  and  experimental  tracking  performance  obtained  with  PI  control, 
optimal  open  loop  control  and  perturbation  control  at  1000  Hz.  (b)  Hysteretic  and  nonlinear  input- 
output  behavior  of  the  magnetic  actuator. 
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•  Department  of  Math.  Sciences,  Montana  State  University,  February  10,  2005. 

•  SPIE’s  12th  Annual  Symposium  on  Smart  Structures  and  Materials,  San  Diego,  CA,  March  9, 
2004. 

•  2005  U.S.  Navy  Workshop  on  Acoustic  Transduction  Materials  and  Devices,  State  College,  PA, 
May  12,  2005. 

•  SIAM  Annual  Meeting,  New  Orleans,  LA,  July  13,  2005. 

•  ICAM  Workshop  on  Mathematics  as  an  Enabling  Science,  Blacksburg,  VA,  October  1,  2005. 

•  ASME  International  Mechanical  Engineering  Congress  and  Exposition,  Orlando,  FL,  November 
9,  2005. 

•  SPIE’s  13th  Annual  Symposium  on  Smart  Structures  and  Materials,  San  Diego,  CA,  February 
28,  2006. 

•  SIAM-SEAS  30th  Annual  Meeting,  Auburn  University,  Auburn,  AL,  March  31,  2006  (Plenary 
Presentation). 

•  CSRI  Workshop  on  Numerical  PDEs  in  the  21st  Century,  Albuquerque,  NM,  April  21,  2006. 

•  Fourth  World  Conference  on  Structural  Control  and  Monitoring,  San  Diego,  CA,  July  12,  2006 
(Keynote  Presentation) . 

•  ICAT  47th  International  Smart  Actuator  Symposium,  Penn  State  University,  State  College, 
PA,  October  4,  2006. 

•  Department  of  Mathematical  Sciences,  Montana  State  University,  Bozeman,  MT,  November  30, 
2006. 

•  45th  IEEE  Conference  on  Decision  and  Control,  San  Diego,  CA,  December  14,  2006. 

•  Department  of  Mechanical  Engineering,  Florida  State  University,  Tallahassee,  FL,  February  1, 
2007. 

•  Idaho  National  Laboratory,  Idaho  Falls,  ID,  March  15,  2007. 

•  Department  of  Mathematics,  University  of  Waterloo,  Waterloo,  Ontario,  April  26,  2007. 

•  Oak  Ridge  National  Laboratory,  Oak  Ridge,  TN,  May  1,  2007. 

•  2007  U.S.  Navy  Workshop  on  Acoustic  Transduction  Materials  and  Devices,  Penn  State  Uni¬ 
versity,  State  College,  PA,  May  15,  2007. 

•  SIAM  Conference  on  Control  and  Its  Applications,  San  Francisco,  CA,  June  29,  2007. 

Transitions 

1.  Magneto strictive  Actuators  -  Etrema:  The  nonlinear  magnetic  constitutive  models  and  pertur¬ 
bation  control  techniques  reported  here  are  being  extended,  in  collaboration  with  scientists  at 
Etrema  Products,  to  optimize  the  performance  of  Active  Machining  Systems  (AMS)  used  to 
mill  products  such  as  piston  heads.  The  goal  is  to  significantly  increase  milling  speeds  (e.g.,  up 
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to  20,000  cycles  per  second)  while  maintaining  micron-level  accuracy.  The  Etrerna  and  AFOSR 
investigations  are  mutually  complementary  and  each  provides  substantial  technology  transfer 
to  the  other.  Point  of  Contact:  Rick  Zrostlik,  Etrerna  Products  Inc.,  Ames,  IA,  515-296-8030. 

2.  SMA  Chevrons  and  Torque  Tubes  -  Boeing:  The  1-D  SMA  models  developed  through  the 
AFOSR  program  are  being  extended  in  collaboration  with  Boeing  scientists  to  2-D  and  3-D 
geometries  inherent  to  chevrons  used  to  reduce  jet  noise  and  decrease  drag  with  potential 
application  to  improved  inlet  channel  design.  Similar  models  are  being  considered  by  Boeing 
as  optimization  tools  for  the  design  of  SMA  torque  tubes  to  change  the  camber  of  rotorcraft 
blades.  In  both  cases,  models  and  control  designs  will  be  validated  using  data  from  Boeing 
experiments  and  flight  tests.  Point  of  contact:  James  Mabe,  Boeing  Phantom  Works,  Seattle, 
WA,  206-655-0091. 

3.  PZT  Unimorphs  -  Boeing:  Nonlinear  structural  models  developed  through  AFOSR  support  are 
being  considered  at  Boeing  for  characterizing  the  hysteretic  and  nonlinear  behavior  of  PZT- 
based  unimorphs  under  investigation  for  flow  control  and  improved  flight  capabilities.  The 
second  phase  of  the  investigation  will  focus  on  model-based  control  design  and  implementation 
of  the  unimorphs.  This  can  potentially  impact  a  broad  range  of  flow  control  problems  of  interest 
to  the  Air  Force  and  Boeing.  Point  of  contact:  James  Mabe,  Boeing  Phantom  Works,  Seattle, 
WA,  206-655-0091. 

4.  SMA  Thin  Films  and  MEMs  -  Sandia:  Shape  memory  alloy  models  developed  through  AFOSR 
support  are  being  investigated  at  Sandia  National  Laboratories  for  characterization  and  control 
design  in  applications  employing  shape  memory  films  and  MEMs.  The  potential  benefit  to  the 
Air  Force  mission  is  significant  since  SMA  films  and  MEMs  retain  the  high  strain  properties 
of  bulk  SMA  but  have  the  potential  for  operating  at  significantly  higher  frequencies.  Point  of 
contact:  James  Redmond,  Sandia  National  Laboratories,  Albuquerque,  NM,  505-844-3136. 

5.  Nanopositioning  -  Asylum  Research:  The  models  quantifying  constitutive  nonlinearities,  hys¬ 
teresis,  thermal  effects,  and  frequency  effects  in  piezoceramic  materials  will  be  employed  in 
conjunction  with  model-based  control  laws  to  improve  the  resolution  and  efficiency  of  nanopo¬ 
sitioners  including  high  speed  scanning  for  atomic  force  microscopy.  Point  of  contact:  Jason 
Cleveland,  Asylum  Research,  Santa  Barbara,  CA,  805-692-2800. 

6.  SMA  Polymers  -  AFRL:  We  are  initiating  investigations  to  extend  the  shape  memory  alloy 
models  developed  through  AFOSR  support  to  shape  memory  resins  and  polymers  presently 
being  considered  by  AFRL  researchers  for  deployment  of  gossamer  structures.  The  SMA  poly¬ 
mers  provide  the  work  density  requirements  necessary  to  deploy  aerospace  structures  such  as 
radar  antennas  and  large  optical  mirrors  but  exhibit  viscoelastic  and  thermodynamic  attributes 
that  are  not  accommodated  by  present  models  or  control  techniques.  Point  of  contact:  Thomas 
Murphey,  Kirtland  AFB,  Albuquerque,  NM,  505-846-9969. 

Publications 

Refereed  publications  resulting  from  work  supported  by  this  grant  are  listed  as  references. 
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